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■^j" , Abstract 

We analyze the proton electromagnetic form factor ratio R(Q 2 ) — 
\»Q ' QF2(Q 2 )/ Fi(Q 2 ) as a function of momentum transfer Q 2 within per- 

■ turbative QCD. We find that the prediction for R(Q 2 ) at large mo- 

mentum transfer Q depends on the exclusive quark wave functions, 
which are unknown. For a wide range of wave functions we find that 
' QFijF\ ~ const at large momentum transfer, in agreement with re- 

C^i cent JLAB data. 

The recent JLAB measurement of the Proton Electromagnetic Form Fac- 
tor ratio [UI21IS1 shows the puzzling behavior R{Q 2 ) = QF 2 (Q 2 )/ Fi(Q 2 ) ~ 
const, at large Q. The form factors F\ and F 2 are defined by the standard 
relation, 
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<p>, s'\J»\p, s >= N(P',s') ^ Fl (Q 2 ) + ia^q v ^p-j N(P,s) (1) 

where q = p' — p, Q 2 = —q 2 , is the electromagnetic current of the proton 
and s and s' refer to the spins of the initial and final proton. F\ and F 2 are 
refered to as the Dirac and Pauli form factors respectively. 

The amplitude iN(p\ s')F2(Q 2 )a )lu q u N(p, s) represents the amplitude 
for chirality of the proton to flip under momentum transfer Q. This flip can 
occur due to flip in the chirality of the proton's constituents. However at 
large momentum transfer Q >> m q , where m q ~ few MeV is the mass of the 
quarks, the amplitude for a quark chirality flip is negligible. Alternatively 
the chirality of the proton can flip due to quark orbital angular momentum 
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[U 03 EE El • General principles allow the existence of substantial quark or- 
bital angular momentum (OAM), but a body of belief accumulated from 
the non-relativistic quark model (NRQM) has not favored this possibility. 
Note that the 3-dimensional OAM of non-relativistic physics has awkward 
Lorentz transformation properties, and we are concerned with OAM of an 
SO (2) subgroup of rotations preserving particle momenta. Here we analyze 
the non-zero quark OAM contribution to Fi in order to determine pQCD 
predictions for the scaling behavior of the ratio R(Q). 

The large Q 2 behavior of form factors is often discussed within the 
Brodsky-LepagejS] factorization scheme. In this scheme the dominant con- 
tributions to amplitudes is assumed to come from the asymptotically short- 
distance (asd) region b ~ 1/Q, where b is the transverse separation between 
quarks. Non-zero OAM is excluded in the first step, and should not con- 
tribute at large Q 2 . The asd formalism is tested by the hadron helicity con- 
servation (HHC) rule [HI A^+A^ = Xc+^D for the reaction A+B — > C+D, 
where A, is the helicity of the particle 'i'. Failures of the asd approach to 
correctly predict experimental results are well known. These include many 
observed violations of the helicity conservation rule and the observation of 
oscillations in fixed-angle proton-proton elastic scattering 10 . A success of 
the formalism is the correct prediction of the Q 2 scaling behavior of many ex- 
clusive processes. Yet scaling can also be obtained from the quark counting 
model of Brodsky and Farrar and Matveev et al , without assuming the 
details of the asd formalism |S]. Thus it is not possible to conclude whether 
or not quarks exist with non-zero OAM by appealing to asd models. 

The elastic scattering of a proton from a virtual photon is shown schemat- 
ically in Figure 1. The three-quark contribution to the proton form factor 
can be written in pQCD as 

N(P',s') ^ Fl (Q 2 )+ia^q v ^^i N(P,s) = J (dk T )(dx)(dk' T )(dx') 

y a /^y(P^k / Ti ,sOr^,^Y a ^ 7 (g,k Ti ,k / T i)Y Q ^ 7 (P,k T i,s) . (2) 

We have factored the amplitude into products of a hard scattering kernel 
and soft initial and final state wave functions Y and Y' respectively. The 
argument k^i of the initial wave function and the hard scattering refers to 
the momenta of the three quarks, kri, kx2 and kr3 with similar definition 
of the argument k Ti of the final state wave function. We use the "brick- 
wall" frame for our calculations. The integration measures are given by, 
(dkx) = d 2 kxid 2 kT2d 2 A>r3<5 2 (kxi + kx2 + kr3), (dx) = dx\dx2dx^,5(xi + 
X2 + X3). Note that Eq. |21 postpones any assumptions about the dominance 
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Figure 1: Schematic illustration of the elastic scattering of a proton from a 
virtual photon, p(P) + 7*(<?) — ► p(P')- 

of any particular integration region of kx or x. If care is not taken, limit 
interchange errors leading to the asd results can result. 

To extract the contribution due to quark OAM, it is convenient to work 
directly with the coordinates b; conjugate to the transverse quark momenta 
kxi- We choose coordinates so that the third quark (down) lies at the origin 
[T2*] . i.e. b3 = 0. The wave function Y a p y can be decomposed as a sum of 
terms [12]: 



We are concerned with terms which lead by power counting in large mo- 
menta. Under a Lorentz transformation along the momentum axes the vari- 
ables b are invariant. We therefore keep leading wave functions whether or 
not a power of b occurs. Each power of 6t —* b x ± ib y can be further decom- 
posed into combinations of quark OAM. (At this point, asd would reject 
powers of b and OAM ^ 0.) The first few operators C are given by 
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(3) 



C, 



l M = (fCU( l5 N) 
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<Ui = KF^cu^n^ (4) 

Here C is the charge-conjugation matrix. Note the operator C 4 which de- 
pends on b\. Here bi are four vectors with transverse components equal to bj 
and all other components zero. Similar operators exist for other transverse 
coordinates. 

The Dirac and Pauli form factors can now be obtained from the Eq. |21 
F\ can be schematically written as 

Fx{Q) = [ (db)(db r )(dx)(dx')<<p*(b, x, x', Q)ip{b, x'). (5) 



Here x, x', Q) is the Fourier transform of the hard scattering 

kernel after projecting out the Dirac bilinear covariant N'j^N, and ip is a 
linear combination of the wave functions V, A and T defined in Eq. |21 A 
detailed calculation of F\ at leading order in perturbation series, neglecting 
the transverse momenta in the Dirac propagators, is given in Ref. 12 . 

We focus on the b integrations to obtain F2 in the limit of zero quark 
masses. Contributing are wave functions such as X (Eq. El), as well as the 
transverse momenta in the Dirac propagators. The scaling behaviour can be 
determined by considering the leading order hard scattering kernel. Since 
we are working in the impact parameter space, we need to take the Fourier 
transform of this kernel. In this case the transverse momentum factors 
such as (&Tl)i> which occur in the Dirac propagators, turn into derivatives 
id/dbu. The remaining Fourier transform has a form similar to what is 
obtained for the form factor F±. The dependence on the impact parameter 



in this kernel arises through the Bessel functions, such as -K"o(y x ix±Qbi2), 

where 612 = |bi — t>2| • Taking the derivative of this kernel with respect to 
bu gives factors of the form, 

(bi-b 

bi2 



Qjx 1 x' 1 K l (Jx 1 x' 1 Q 2 b 2 12 



Besides this transverse separation dependence, an additional power of b 
arises from the operator C 4 in the wave function. 

The scaling of F2(Q 2 ), and the form factor ratio, hinges on scaling of 
the effective transverse separation b at large Q. Counting powers of Q, 
including one for the prefactor q^a^v, we find that if 6 ~ 1/Q at large Q, 
then F2/F1 ~ l/Q 2 in this limit. This result has recently been confirmed 
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in Ref. ^3] which adopts the asd formalism from the start. Alternatively 
if b ~ constant then QF2/F1 ~ const and F2 ~ 1/Q 5 in the same limit. 
Such scaling was predicted in Ref. [5] and is also seen in a relativistic quark 
model calculation |14j . 

We turn to how the dominant region of b actually scales with Q. The 
simpler case of the pion form factor [To] is instructive. The one gluon ex- 
change kernel in this case can be written as 



H(b,x, x', Q) = 87r 2 C F a s K (^xx'Q 2 b 2 ). (6) 

where a s is the strong coupling and Cp is the color factor. At large Q the 
dominant contribution is obtained from the region 

Vxx'Qb < 1. (7) 

In order to reproduce the asd assumptions, average values of x ~ 1/2 for pion 
and x ~ 1/3 for proton can be assigned, converting V xx'Qb < 1 — > Qb < 1. 
Such estimates are not the same thing as actually computing the result! 
Indeed one cannot rule out the possibility that in the limit of very large 
Q, b ~ constant and V 1 xx 1 ~ 1/Q, or any intermediate combinations of 
integration regions, given the complexity of the problem. Elsewhere [7] we 
have shown that it is not possible, in principle, to determine the scaling 
of a power of b without specifying the interplay of x and 6-dependence in 
the wave functions. Indeed the "rules of pQCD" require one to objectively 
treat wave functions as unknown non-perturbative objects, to be determined 
from data. Constructing wave functions to make the asd results come out 
is possible, but we find such logic circular. 

Suppose one uses a Gaussian model for the wave function x, x') = 
b A e~ b2 /( 2a2 ^<fi(x, x') to represent the wave functions cutting off large b. The 
factor of b A is the phase-space to find A quarks close together from naive 
quark-counting. We probe effects of quark OAM by calculating the moment 
< b(Q) > T defined by 

J dxdx'd 2 bbJ r 7r (Q, x, x' , b) . , 

< HQ) >.= . (8) 

Then the pion form factor F 1T {Q 2 ) is given by 

FAQ 2 ) = J dxdx'd 2 bF 7T (Q, x, x', b) (9) 

and J- n = H<&. We change variables for the longitudinal fractions to £ = 
V xx', ( = x/x' and parameterize </>(£, C) ~ £ r+1 (l — <!;) r+1< MC)> as £ — > 1: 
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the £ dependence can be left unspecified. By substituting this into Eq. 
|H1 when r < A, the dominant contribution to the numerator in Eq. |H] is 
obtained from x ~ 1/Q, and the moment < b(Q) > 7r ~' const. Due to the 
inherent unknowns of the wave function's dependence on both x and b, it is 
not possible to establish suppression of the contributions of quark OAM in 
pQCD in a model-independent way. 

Similar results are obtained for the proton. In Fig. Owe plot one moment 
of the transverse separation bi for the proton using the COZ model for 
the x-dependence of wave functions [16] . both including and not including 
the Sudakov form factor [15] ■ It is clear that over a very large range the 
moment has a very weak dependence on Q. Hence the form factor ratio 
R(Q) = QF<ijF\ scales like a constant in this region, in agreement with the 
recent JLAB experimental result. Asymptotic analysis shows that for a 
wide range of x-dependence in wave functions, the form factor ratio R(Q) 
scales like a constant, modulo logarithms. 
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Figure 2: The moment < 62 > of the proton form factor kernel with and 
without including the Sudakov form factor. The COZ model is used for the 
x-dependence of wave functions. 
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